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QUANTUM HOMODYNE TOMOGRAPHY AS AN INFORMATIONALLY 
COMPLETE POSITIVE OPERATOR VALUED MEASURE 

PAOLO ALBINI, ERNESTO DE VITO, AND ALESSANDRO TOIGO 



Abstract. We define a positive operator valued measure E on [0, 2tt] x K describing the measurement 
of randomly sampled quadratures in quantum homodyne tomography, and we study its probabilistic 
properties. Moreover, we give a mathematical analysis of the relation between the description of a state 
in terms of E and the description provided by its Wigner transform. 



Quantum homodyne tomography; positive operator valued measures; Wigner function. 
CN ' 1. Introduction 



Quantum homodyne tomography [H [21 [3l [4] allows to determine the state of a single mode radiation 
field by repeated measurements of the quadrature observables Xg, the phases 9 being chosen randomly 
in T = [0, 2tt]. This can be seen as a consequence of the fact [5; that, for a large class of observables O, 
'^ ' there exists an associated function /o : T x R ^ R such that 



(1) tv[Op] 



fo{0,x)duPg{x) 



de 

2^' 



where Vg is the probability distribution of outcomes obtained for the quadrature Xg measured on the 

^ ■ state p and d0/27r is the uniform probability distribution on T. Actual reconstrunction schemes are 

f — I strictly related to a statistical interpetation of formulas of this kind. Indeed, quantum tomography 

m . experiments output a n-uple {(8i,Xi)}"^^ of pairs in T x M, each one of which represents the outcome 

Xi of a measurement of the quadrature observable corresponding to the randomly picked phase 0^. If 

one assumes such pairs to be samples from a random variable on T x R distributed accordingly to a 

r^ I probability measure /i'' such that 

o: (2) dp''{e,x) = di.^g{x)^, 

one can use the experimental outcomes to estimate integrals such as (HI) (see Ref. [2j and references 
therein), for example by replacing d/i'' with its empirical estimate — X]"=i ^iQi,Xi)- 
}J] ' The above reconstruction formula, although very popular, is not the only scheme used for tomograph- 

5^ \ ical state estimation: other ones are known which don't rely on it [2]. The hypotesis that experimental 

results are distributed accordingly to ([2]) lies however under both every proposed reconstruction algo- 
rithm and its statistical analysis [U [2l [S] [6l [7] . Actually, although given for granted in the cited literature, 
well-definiteness of a joint probability distribution such as /i'' in ([21) is a priori not trivial. In the first part 
of our paper we prove well defineteness of fj.^ by showing there exists a positive operator valued measure 
(POVM) ^ on T X R such that 

pP{Z) = tr[E{Z)p] for any Borel subset Z of T x M. 

According to the physical meaning of /i'', the POVM E is the generalized observable associated with the 
quantum homodyne tomography experimental setup. In particular, we show that /i'' has density p''{9^ x) 
with respect to the Lebesgue measure on T x R, its support is always an unbounded set, and the mapping 
p i-^ pP \s injective (i. e. E is informationally complete). The intertwining property Xg — e Xe^^^^ , 
where N is the number operator and X is the position operator, turns out to be crucial for the definition 
of E (or, equivalently, for the definition of /i''). We remark that the introduction of a POVM for the 
homodyne tomography measurement process is already present in physical literature (see section 2.3.2 
in Ref. [2]), but it is grounded on a formal construction. We provide here an alternative, rigorous 
formulation. 

In their seminal paper on quantum homodyne tomography [4], Vogel and Risken argued that the the 
Radon transform of W{p), where W{p) is the Wigner function associated to p, is precisely the proba- 
bility density function p^ generated by the homodyne tomography measurement, so that the following 
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commutative diagram holds: 

States on H. 




Radon transform 

Wigner functions on M ^ probability densities on T x M 

The suggested estimation procedure, applied also in the first homodyne tomography experiments, is then 
based on the inversion of the Radon transform by means of classical techniques in medical tomography. 
However, the derivation of this fact is once again rather formal, and never given a rigorous basis in the 
literature on the subject: in the second part of the paper, we will thus address the problems that arise 
in looking at such formulation of quantum tomography from a rigorous point of view. First of all, we 
will recall that in order for the Radon transform to be well-defined, we need the Wigner function W{p) 
to be integrable on M^. Then we will show that the support of W{p) can never be bounded. This is a 
potential problem, since the estimation techniques used in classical tomography are explicitly devised for 
compactly supported objects. One can however by-pass the problem and still give an inverse for the Radon 
transform if he assumes that the Wigner function under observation is a Schwartz function on M^ . This is 
precisely what happens in most homodyne tomography experiments, where the states under observation 
are linear combinations of coherent or number states. In Section [4] we show that this assumption on the 
Wigner function is equivalent to suppose that p has a kernel which is a Schwartz function on R^ (since 
p is an Hilbert-Schmidt operator, p is an integral operator whose kernel is a function on R^). Under this 
assumption on p we prove that the Radon transform of W{p) is p^ and the inversion formula holds true. 

2. Preliminaries and notations 

In this section, we will introduce the notations and give a very brief description of the mathematical 
structure of quantum homodyne tomography. 

2.1. Notations. Let 7i be a complex, separable Hilbert space with norm ||-|| and scalar product (• , •) 
linear in the second entry. Denote by C{Ti.) the Banach space of the bounded operators on Ti with uniform 
norm ||-||£. Let Xi(7i) be the Banach space of the trace class operators on Ti. with trace class norm ||-|j-^, 

and let 5(7i) be the convex subset of positive trace one elements in Xi(7i). Finally, let 12{'H) be the 

1 /? 
Hilbert-Schmidt operators on Ti, with norm WAW^ — [tr [yl*yl]] ' . We recall that the elements of 5(7Y) 

are the states of the quantum system whose associated Hilbert space is TL. 

Suppose 51 is a Hausdorff locally compact second countable topological space. Let B{^) be the Borel 

cr-algebra of fi. We recall the following definition of positive operator valued measure. 

definition 1. A positive operator valued measure (POVM) on D, with values in 7i is a map E : B{^) — > 
C{Ti.) such that 

(i) E{A) > for aU A e B{n); 
(ii) Bin) = /; 
(iii) if {Ai}i^i is a denumerable sequence of pairwise disjoint sets in B{Q), then 

E{ua^) = y^.EiA,), 
^ — ^i 

where the sum converges in the weak (or, equivalently, ultraweak or strong) topology of CiTL). 

E is a. projection valued measure (PVM) if E{AY = E{A) for all A e B{Vl). 

If £; is a POVM and T e Ii(H), we define 

p%{A)^tY[E{A)T] VAe6(f^). 

Then, prj^ is a bounded complex measure on 57. If p G S{Ti), /i^ is actually a probability measure on 
ri, and /i^(A) is the probability of obtaining a result in A when performing a measurement of E on the 
state p. 

2.2. The mathematics of quantum homodyne tomography. The physical system of quantum 
homodyne tomography is a single radiation mode of the electromagnetic field. The associated Hilbert 
space is H = L'^ (M). Let 



^ = -^ p{x) e 2 I p is a polinomial, > 



QUANTUM HOMODYNE TOMOGRAPHY AS AN INFORMATIONALLY COMPLETE POVM 3 

which is a dense subspace of 7i. As usual, we denote by X and P the position and momentum operators, 
respectively. Their action on A is explicitly given by 

{Xf){x) = xf{x) and {Pf){x) = ~i^{x). 

ax 

Letting T = [0, 27r], for any ^T the corresponding quadrature is the self-adjoint operator Xg on L? (M), 

whose action on A is 

Xg = coseX + sm9P. 

li x,y €z M, and x — r cos 6, y = r sin 9, we have 



[e'"-^V] (2) = \e''^^^+yP'>f] (z) = e<'^+"")/(z + y) 



for all/ e L2(R). 

We denote by Hg the PVM on R associated to Xg by spectral theorem. In particular, 11(^4) := no(^) 
is just multiplication in L^ (M) by the characteristic function 1a of A, while niL(j4) — !F*Il{A)T, where 
J-' is the Fourier transform 

(3) Tf=^ [ e-^^yf{y)dy f e L' L^R). 



The number operator is the essentially self-adjoint operator N whose action on A is 

N^ -{X^ + P^ -I) . 

For all 6* e T, we let V{9) — e'^^. Since the spectrum of N is N, the map 9 -^ V{9) is a unitary 
continuous representation of T acting on L^ (R) , where we regard T as a topological abelian group with 
addition modulo 27r. The number representation V intertwines the quadratures Xg, in the sense that 

Xg = V{9)XV{9)* 

for all 9 eT, and 

ng{A) = v{9)n{A)v{9y 

ior all 9 e T and A & B{R). 

Finally, given p G S{'H) and 9 G T, we denote by Vg the probability distribution on M of the outcomes 
of the quadrature Xg measured on the state p, namely 

(4) i/^(A) = tr [pUg{A)] = tr [pV {9)U{A)V (9)*] VA e B{R). 

3. Main results 

In this section, we will describe explicitly the POVM which intervenes in homodyne tomography and 
the associated probability distributions on states. 

The first result studies some properties of the family of probability measures Vg defined by ([5]). In its 
proof and in the statement of some of the following results, we will make use of the concept of section 
through some G T of a Borel set S G S(T x R), defined as follows: 

B^ = {.TeR|(6l,a;) G B}. 

Proposition 1. Given p E S{T-i) and 9 E T 

(i) the probability measure Vg has density pg G L^ (R) with respect to the Lebesgue measure on R; 
(ii) the map 9 i— *■ Vg{B ) is measurable for any B E B{T xR). 

Proof, li A E B{R) has zero Lebesgue measure, then lV{A)f = 1^/ = for aU f E L'^ (R). Therefore, 
vPg{A) = tr [pV{9)\i{A)V{9)*] = 0. Thus, the first claim follows. 
If {cnlnsN IS a Hilbcrt basis of Ji, then 

tv[pV{9)Tl{B')V{9Y] = ^ (e„ , y(0) Vy(0)n(i3«)e„) 

n 

= im (en - V{9ypV{9)e„,) {era , n(B'')e„) . 



Since the map 9 i— s- (e„ , V{9y pV{9)em) is continuous and the map 9 i-^ (cm , n(i?^)e„) = f Ib{&, x)en{x)em{x)dx 
is measurable by Fubini theorem, measurability of 6* i-^ tr [pV {9)Il{B^ )V (9)*'^ follows. D 

Next theorem shows the existence of a POVM associated to quantum homodyne tomography. This 
theorem should be compared with the formal derivation of E given in Ref. ^ (see eq. (2.34) therein). 
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Theorem 1. There exists a unique positive operator valued measure E on T x M. acting in L^ (M) such 
that 

(5) tr[pE{B)]= j^uP,{B%^ 

for all p e S{n) and B e B{T x R). 

Proof. Eq. (g]) suggests to define the POVM as 

E{B)^ f V{d)n{B'>)V{9)*^. 

Jt 27r 

To prove that the above definition is correct, we first show that the map i-^ V {9)I1{B^ )V {6)* is ^- 
ultraweakly integrable for aU Borel subsets S of T x E, and then we prove that B i-^ E{B) is a POVM. 
Now, given B e B{T x R) and p e S{n), the map ^ tr [pV {e)n[B%V {6)*] is measurable by the 
previous proposition, and 

|tr [pV{e)i].{B'')v{e)*] I < iIpIIi ||n(B^)||^ < i ve* e t. 

Therefore, it is |^-integrable. This shows that 9 ^ V {e)Il{B^)V {9)* is |^-ultraweakly integrable. 

Suppose T e Xi(H). Then T = ELo^'^fe' ^i^h T, > and liTolli + {t^W^ = ||Ti||i + HTslli < WUv 
Setting pk = Tkl \\Tk\\-y (with 0/0 = 0), we see that 



tr [TV{9)I{{B'^)V{9)*] — 



A9 ^ 



J2\\n\\J \tT[p,v{9MB')v{9r]\- 



< 

fc=0 
3 

< 5]||Tfc||,<2||r||,. 

fc=0 

This shows the existence of E{B) e C{H). Clearly, E{B) > 0, and E{T x M) = /. 

If {i?„}„gN is a monotone increasing family of elements in B{T x R), with Bn t -B, then, for all 9, 

tr [pV{9)UiB'jVi9r] = ^,^(i?^) T ^'eiB') - tr [pF(0)n(i3«)F(0)*] . 

By dominated convergence theorem 

r\f) f c\f) 

tr[pV{9)n{Bt)V{9r] - t / tr[pV{9)n{B<^)V{9r] -, 

and ultraweak cr-additivity of E follows. D 

We let pP — tr[_E(-)p] be the probability distribution on T x R associated to a measurement of E 
performed on the state p. By definition (O it follows that 

)^ 

/T 27r 

as wanted. The following theorem gives some properties of p''. 

Theorem 2. Let p e S{n). 

(i) The measure p^ has density with respect to ^dx. We denote such density hy p^ . 
(ii) For ^-almost all 9, p''{9, x) — Pg{x) for Ax-almost all x. 

(iii) The marginal probability distribution induced by p^ on T is the Haar measure ^, and the condi- 
tional probability distribution induced by pP on R is Vg for ^'-O'^^^ost all 9. 

Proof (i) If B e B{T x R) is a || dx-nuU set, then B^ is dx-nuU for ff -almost all 9 by Fubini 

theorem, so, for such 6''s, i/g{B^) = 0. Therefore, p''{B) = by ([6]), thus showing that pP has 
density with respect to ^ dx. 
(u) If Z e 6(T), A e S(R), we have 



(6) p^iB)^ I <{B')\ 

Jt 






This holds for all Z, implying that there exists a ^-nuU set Na G B{T) such that p''{9,-) is 
dx-integrable with 

[ pPi0,x)dx^,yP{A) 
Ja 
for all e^ Na- 
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Let {An}neN be a sequence in S(R) with the foUowing property: if /ii, /i2 are positive measures 
on M such that /ii(A„) = /i2(A„) for all ti, then iii = /i2 (such sequence exists since M is second 
countable by Theorems C §5 and A §13 in Ref. [SI). Let A^ = U„iVA„- Then N is ff-nuU, and, 
ii ^ N, pP{9, •) is integrable with 



p^{6,x)dx— / pg{x)dx Vn. 

This implies that, if 6* ^ TV, p''(0, a;) = Pg{x) for dx-almost all x. 
(iii) This is just dH). 

D 

Remark 1. As a consequence of item (iii) in the above proposition, a well known result on conditional 
probability distribution ensures that, if </) is a /i^'-integrable function, then 0(0, •) is i^g-integrable for 
l^-almost all 9, the map 6 i-^ J^(j){9,x)di^g{x) is ^-integrable, and 

/ 0(0, x)dfiP{0, x)^ j\j 0(0, x)dvP{x)\ ^. 

By Theorem [21 E is the POVM associated to the measurement of a quadrature Xq chosen randomly 
from T with uniform probability ^. 

The next corollary shows that the probability distribution /i'' can not have compact support for any 
pizSi'H). 

Corollary 1. For all R > 0, we have 

[ [ pP{0,x)dx^>O. 

JtJ\x\>r ^tt 

Proof. With AR = {xeR\ \x\ > R}, we have 

/ / pP{0, x)dx^ = ^,f{T xAn)^ f tr [pV i0)niAR)V (0)*] ^ 

= tv[p'U{AR)], 

with p' = Jrj. V{0)* pV{0)^. p' is a trace one positive operator. Since it commutes with the representation 
V of T, it is diagonal in the number basis {e„}„gN of L^ (M). Since (e„ , n(Afl)e„) > for all n, the 
claim follows. D 

As a consequence, the map p i-^ p^ from 5(7i) to the set P(T x M) of probability densities in L^ (T x K) 
is not surjective. The next corollary shows that it is actually injective, i. e. the POVM E is informationally 
complete [9]. 

Corollary 2. If p,a- E S{'H) and p ^ a, then pP "^ p" ■ 

Proof. If p, cr G S{T-C), then pP = p'^ if and only if pP — p"' (in L^ (T x R)), which amounts to say that 
Pg — Pe (i^ ^^ (^)) for ^-almost all 0. This is in turn equivalent to Vg — Vg for |^-almost all 0. For 
r e M and € T, we have by spectral theorem 



*™di^^(a;)= / e'™tr[pne(da;)] ^trLpe'"--^"] = V27r[y(p)](rcos0,rsin0), 
where 



[V{p)]{x,y)^^iv 



p^t{xX+yP) 

Since the map V : Ii(H) — > C(R^) is injective (see for example Ref. tlQi), injectivity of the map p i-^ pP 
follows. D 

4. The Radon transform of the Wigner function and Radon reconstruction formula 

In the previous section, by means of the POVM E defined in Theorem [1] we estabilished a convex 
injective correspondence p ^^ pP between states and the set of probability densities on T x M. However, 
no explicit formula relating p to the function pP was given, due to the fact that, if p does not have a 
simple expression in terms of the number basis, the expression tr [T^(0)*py(0)n(^)] can not be explicitly 
computed. 

In this section, we will show that, if the state p is sufficiently regular, pP can indeed be evaluated, being 
in fact the Radon transform of the Wigner function W{p) of p. This is a very well known fact in quantum 
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tomography, going back to the seminal paper of Vogel and Risken j4j. However, no attention has never 
been paid in the hterature to the fact that performing the Radon transform of W{p) makes sense only 
for a restricted class of states, namely for those p G S{7i) such that W{p) G L^ (R^). This constraint 
becomes even more stringent when one considers the inverse formula reconstructing p (or, better, W{p)) 
from its associated probability density p'^. We will see that, in order to derive mathematically consistent 
formulas both for pP and the reconstruction of W{p), one needs to assume that the state belongs to the 
set of Schwartz functions on M^. This seems a rather strong limitation, as the very natural attempt to 
extend the Radon transform and Radon reconstruction to the whole set 5(7i) by means of distribution 
theory fails in the quantum context (Remark [2]). Our main reference to the results below is Ref. |11) . 
If T G Xi(7i), we introduce the bounded continuous function V{T) on E^, given by 



(7) [V{T)]{x,y) = ^tT 



rp^i(xX+yP) 



It is well known (see for example Ref. (TU]) that V{T) e L^ (K^), and V uniquely extends to a unitary 
operator V : 2'2(7i) — > L^ (]R^). The Wigner transform of ^ G 2^2 (7i) is just (up to a constant) the 
Fourier transform of V{A), i. e. 

(8) W{A)^^T,V{A), 

y Zn 

where T2 = J^ (S> :F on L'^ (R^) = ^2 (^^-^ ^ ip. (]j)^ ^-^j^ j: defined in ©. 

If / G L^ (K^), the Radon transform of / is the complex function Rf G L^ (T x M) given by 

/+00 
f{r cose~t sin 6*, r sin + t cos e)dt 
-00 

^ dr-almost everywhere. 
We have the following fact. 

Proposition 2. IfW{p) G L^ (M^), then 

(10) [RW{p)]{e,r)=pP{9,r) 

for ^dr-almost all {0,r). 

Proof. Let 7 : T X K — >M?he the map 

7(0, r) — (r cos 9, r sin 9) . 

We have 

[V{p) o 7](0, r) = ^=tr [pe-^«] = ^= / e-*pg(t)dt = [T^Vel (0 

V Ztt y/ Ztt J -00 

by spectral theorem. On the other hand, 

1 /•+00 (•+00 
[T^^W{p)oj]i9,r)^— / e^(-'-™^^+^'^-"«)[W^(p)](x,2/)dxdy 

^J_f f e"'^(™" '^""'^ *+'''" '^'''"^)[W^(/9)](t cos (/),t sin 0)|i|dt^ 
27r Jo J -00 27r 

^TT Jo J -00 ZTT 

^^fl^"°e^'''^°''^[W{p)]{tcos{4> + 9),t8m{cl, + 9))\t\dt^ 



1 n+oc ^+00 

— / / e*™[VF(/9)](a;cos0- ysin^, ycos0 + xsin0)dxdj/ 



2tt 
1 



00 
+ OC 



2^ , e^^'^[RWip)]{9,x)dx 



' —OQ 

1 



T-'[RWip)]i9,-)]ir). 



V27r 
By injectivity of Fourier transform, the claim then follows by comparison. D 

Corollary 3. The support ofW{p) is an unbounded subset o/R^ for all p G S{H). 
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Proof. Suppose by contradiction that W{p) = almost everywhere outside the disk Dr of radius R in 
M2. Then W{p) G L^ (M^), and so [RW{p)]{e,r) ^pP{e,r) by the above proposition. We have 



|r|>fl 



\RW{p){e,r)\dr^ 



< 



'\Da 



do 

\W{p){rcos9 — tsm9,rsin9 + tcos6)\drdt — 

27r 



'■\Dr 



do 

\W{p){r,t)\drdt— 



0, 



which contradicts CoroUary[TJ 



D 



The first formal derivation of PH)) is contained in Ref. Pj, without the assumption W{p) G L^ (M^). 
We stress that if W{p) ^ L^ (K^), then pi7| does not make sense, and the only possible definition of p'' 
is by means of item 1 in Theorem [2j 

If we denote by S^{T-i) the subset of states p £ S{Ti) such that W{p) £ L^ (K^), then we have 
estabilished the following diagram 



^ P{T X M) 



Now we turn to the problem of reconstructing W{p) given pf. If W{p) € S{R'^), the space of Schwartz 
functions on R^, Radon inversion formula is applicable, and we can obtain W{p) from pP in a rather 
explicit way. Before stating Radon inversion theorem, according to Ref. [llj we need to introduce the set 
ShQ!'^) of functions : T x M — » C such that 

(i) 0e C°°(TxM); 




< oo; 

7" 



(ii) sup,^, (i + |^|fc)^^0(0,^) 
(iii) 0(6*, r) = 0(27r - 0, -r) for aU 
(iv) for each k E N, J_ <j){9, r)r'^dr is a homogeneous polynomial in sinO, cos 9 of degree k. 

It is shown in Ref. [H] that Rf £ Sh{F^) if / G S'(M^), and the map R : SiW^) — > Sh{P^) is one-to-one 
and onto. Thus, in our case W{p) G S{M.'^) is equivalent p'^ £ Sh(P'^) by Proposition [2l 

The next theorem is a restatement of Theorem 3.6 in Ref. [IT] (see also Ref. 4 for a formal derivation 
of pTjl ). We stress that the hypothesis W{p) £ 5(M^) (or, equivalently, p'' £ ShO^"^)) is needed in order 
to give meaning to p2|) and to define the integral in p^ . 



Theorem 3. Suppose W{p) £ S{M.'^). Then 



(11) 
where 

(12) 

and 

(13) 



W{p) 



^R*l^P^ 



ApP{9,r) 



[Ttm]*p''{9,-)]{r)^PY 



+ OC 



1 dpP{9,t) 



r-t 



dt 



dt 



R*f{x,y) 



27T 



f{9,xcos9 + ysm9)— V/ G C°°(T x 

ZTT 



(in (|12p . the Fourier transform of \t\ and the convolution are interpreted in the sense of tempered distri- 
butions, and PV is the Cauchy principal value of the integral). 



We devote the rest of this section to find the subset of states p £ 5(7i) such that W{p) £ 
i. e. to which both Radon transform (jlOp and Radon reconstruction formula (jlip are applicable. 

Each T £ Xi{'H), being a Hilbert- Schmidt operator on L^ (M), is an integral operator, whose kernel 
Kt is in L^ (M^). We have the following fact. 



Proposition 3. Suppose K £ 
trace is 

(14) 



^). Then the integral operator L^ with kernel K is in Tiiji.), and its 
K{x^x)dx. 

-oo 
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Moreover, Lk G 'S{'H) if and only if K is positive semidefinit\} and j_ K{x,x)dx = 1. 
Proof. Let / = (— tt, tt), and let $ : L^ (M) — > L^ (/) be the following unitary operator 

$/(j/) = (l + tan2y)l/2/(tany). 
$ intertwines Lj^ with the integral operator L^ on L^ (/) with kernel 

K{yi,y2) = {l + tav?Vif^K{t&nyx,ia.ny2){l + t&v?V2Y''^ 2/i,2/2 e (-7r,7r). 



Since iiT extends to a C°°-function on / x / by setting iiT = in the frontier of / x /, by Lemma 10.11 



in Ref. [12] L^ is a trace class operator on L^ (/), whose trace is given by 



/7r /'+OC 

-TT >^— OO 



Since L^f = <I)~iL^$, eq. dH]) follows. 

It is easy to check that, if K is positive semidefinite, then the integral operator Lk is positive. 
Conversely, suppose Lk is a positive operator. Fix a Dirac sequence {/„}„gN, and let gn = X]i=i Ci/nS 
where fn^{x) = /„(a; — x,;). We have 

< {gn, LKgn) = ^ c,cj / /„(a;-Xj)i\:(x,j/)/„(y-a;i)da::dy 

- ■ 1 J — oo J — oo 



ij=l •'-■^ •'-'^ 



Af 



n — *oo 



/ ^ CiCjI\ [^Xj , Xi), 



from which positive definiteness of ii' follows. The last claim in the statement is thus clear. D 

We introduce the following linear subspace of Xi (H) 

If in) = {r e x^{n) \ Kt e s{R^)} , 

and define 

If r e Xf (7i), we can explicitly evaluate the trace in ([7]) and the Fourier transform in ^ defining V{T) 
and W{T) respectively. We find 

[V{T)]{x, y) = T^^ [Krit + y/2, t - y/2)] {x) 
[W{T)]{x, y) = -^Tt [Kt{x + t/2, x - t/2)] (y). 



where we denoted by Tt the Fourier transform with respect to the variable t. The second formula proves 
the next proposition. 

Proposition 4. W : if (W) — > SiM?) is a bijection. 

Restricting to states in S^{Ti.), we have thus arrived at the following diagram. 




Remark 2. Unfortunately, one can not use the definition of Radon transform of distributions to extend 
([T0|) to whole L^ (R'^) , or reconstruction formula (fTTj) to a larger set than S^{Tl). In fact, as explained in §5 
of Ref. [11] , the distributional Radon transform can be defined only as a map R : £'{TxM.) — > £' (T x M) , 
£'{T X R) being the set of compactly supported distributions on T x M. Corollary [3] then prevents us 
from giving any distributional sense to (jlOp . Similarly, eq. pip has no distributional analogue, as the 
reconstruction formula T = j^R^[ART] (Theorem 5.5 in Ref. JJJ) again holds only for compactly 
supported distributions T. 



We recall that a function K : M.'^ > C is positive semidefinite if ^^ --^ CiCjK{xj,Xi) > for all Af g N, ci, C2 . . . cjv C 

C and xi , a;2 . . . xj^ C IR. 
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Remark 3. Being able to exhibit an explicit inversion formula of the Radon transform only for Wigner 
functions which are Schwartz class does not imply a failure of quantum tomographical methods in re- 
constructing states with weaker regularity properties, as the associated POVM remains informationally 
complete on the whole of S{Ti.), as we have shown in the first part of this paper. In fact, mainly in order 
to address issues of numerical stability, actual reconstruction methods usually do not involve -^R^K 
directly, but some approximated technique involving regularizations; proofs of consistency are available 
[T] for some of these regularized estimators which holds on the whole of quantum state space. 
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